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ABSTRACT. In this work, we present reproducing kernel method to solve non- 
linear system of differential equations. Several numerical results are given to 
show the efficiency of the method. 
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1. INTRODUCTION 


System of differential equations and dynamical system have been widely used in 
applications such as mathematical model, statistics, biological application, popula- 
tions etc. Consider a system of differential equations 


w(t) = f(t, x(t), y(t)); 


y'(t) = g(t, x(t), y(t); 
subject to 2(0) = 0 = y(0), where f and g are continues in t, x, and y. 


(1.1) 


2. PRELIMINARY 


This section is devoted to some basic definitions, lemmas, and theorems which 
will be used in later. 


Definition 2.1 ( [1]). Let H be a Reproducing kernel Hilbert space of continuous 
real-valued functions f : S —> R. Then there exists k: S x S > R in H such that 
(fC), kC, 8))a = f(s), for all f €H and all se S. 


Definition 2.2 ( [1]). Let W3[0,7] be the space of all functions u : [0,7] > R 
such that u is absolutely continuous on [0,7] and its derivative belongs to L7[0, T]. 


Definition 2.3. Let W2[0,T] be the space of all functions v : [0,7] > R such that 
v’ is absolutely continuous on [0, T] and v” belongs to L?[0,7']. 


Definition 2.4 ( [2]). The space W4[0, T] equipped with the inner product 
(uho = m(O)0(0) + fal (rola 


is a RKHS with reproducing kernel function q(-,-) given by 


(2.1) q(T) = cosh(r +t- T) + cosh(| 7 -t| -T)). 


way! 
2sinh(T) 
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Definition 2.5. [1,2] The space W2[0, T] equipped with the inner product 
(ue) gqocry = WOJO) +(+ fw" (ryu"(r)ar, 
[0,7] 


is a RKHS with reproducing kernel function r(-,-) given by 


(2.2) later 


Definition 2.6. [2] Let W’[0,T] = W2[0,T] 
vector functions fi : [0,7] x [0,7] > R where 
for j = 1,2. 


© W3[0,T] be the space of all real 
h = [hi,h2]" and hy, ho € WE[0,T] 


3. REPRODUCING KERNEL METHOD 


In this section, the general structure of representation of the exact solutions will 
be discussed to problem (4.1). 

We define operators F,G : W3[0,T] + W3[0,7], for each x,y € W3[0,T] such 
that 


| Fe()=2"() 

(3.1) 

Gy) =y"). 

System (4.1) can be written as: 

| Fa(-) = fC #0), 90) 
Gy(-) = 9(,#(), ¥C)). 


Define an operator A: W?[0,7] > WH[0, T], such that systems (3.1) and (3.2) 
can be represented as: 


(3.2) 


=> 


(3.3) AU =w, 

where A = [F,G]" and w! = [zy]. 

Lemma 3.1. The operators F,G : W2[0,T] > W3[0,t] are bounded. 
Proof. Let x €e W3[0,T]. Then 


2 
[Eal [0,T] -(Fe, Pa) waren 


Pena 
T 
=)? f CAO 
T 2 
< (C0)? eo)? f CAO. 


2 
= læ wW2[o,T]' 
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Similar argument can be used to show that the operator G from W2[0, T] into 
W3[0, 7] is bounded. 

Definition 3.2. Let j = 1,2. The space W’[0,7] equipped with the inner product 


X,Y 


? 


he. ta (£1, Y1)w [0,T] + (22,2 lwz [0,7]? 


isa RK HS with reproducing kernel function Qas) = [a(-),¢s(-)]' where q,(-) 


is defined in (2.1) and Rea) = [ri(-), rs(-)]" where r,,(-) is defined in (2.2), 
respectively. 


Lemma 3.3. The operator A: W°[0, T] > W'[0,T7] is bounded linear operator. 


Proof. Let X= [21,22]" belongs to W?[0,7] where x1, 22 € W3[0,T]. Then 


= )2 
(e.g Peers =|X Renan 


=(x’, ee 


=(21 501) [0,7] F (x9, ©) wa to.7) 
a 2 1 2 j " 2 
=(2',(0))2 + f (2'(4))Pat + («4 (0))? + T COL 
T 2 E 2 
<(1(0))? + (21 (0))? + i: (2! (t))at + (22(0))? + (24(0))? + CO 
0 0 
=lz1 lwo r]+lz2lw2to,r] 


HX lwar 


Let {(&,7;)}#j-1 Þe a countable dense subset of D = [0,T] x [0,T] and define 
Wi; = [hip] = NO aces where A* is the adjoint operator of A and Oech is 
the kernel function of W'[0, T]. 


Lemma 3.4. Let {(&,7;)}?°_, be a countable dense set of points in D. Then Vi 


i,j=l 
is a complete system in W7[0,T], and 


Vij =ARG,,s;) 
for all i,j €N. 


Proof. Let Ñ = [x,y]" € W?[0,7] and let g e a 0, for all i,j € N. 
Then ' 


(zY ta ie oT] 


4G Ois(s, ae 
=A (s;,tj) = 0, 


w![0,T] 
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where U (s;,tj) = [x(si),y(t;)]". Since AW@(s;,t;) is vanish in a countable dense 
set of points in D, and it is continuous in W1[0,T] then w = 0. Furthermore, 


Wig =A" Q Céir) 
-(A* Q (éit) C, ‘), Hea; Gas, T] 


-( Qir) J), ARs; Jhe 
=AR¢,,s;)- 


*[0,T] 


The proof is complete. 


In fact, Gram-Schmidt orthonormalization process can be used to drive normal 
orthogonal basis {9 ;}%, of W?[0,T] from {V;;}%, as follows: 


i T 
Oi | Bike, >, aut 
k=1 k=l 
i Ex z 
=> Bir V rk = [8:8] , 
where Big are orthogonalization coefficients of Yp. 


Theorem 3.5. Let {(£i,7;)}7j-1 be a countable dense subset of points of D and 


let X € W?[0, T] be a solution of AX = X is defined in (8.8). Then the solution 
X has the following representation: 


(3.4) RF Y Bux’, 


i=1 k=1 Ww?[0,T] 
= YA AX, Be, w ®; 
i=l k=1 W![0,T] 


3° Y BAX (te) ®; 
k=1 


=> $ bX’ (th) Bi. 
i=1 k=1 
Since X’ (tp) = [x (tk) y'(te)] = [tr £k: Yk), ltr, ks Ye)] , so the proof is com- 


plete. 
4. NUMERICAL RESULTS 


Example [3] A Lotka Volterra model is described by the following system of 
differential equations: 


de = ax — bry 
(4.1) 
dy = -cy + dxy 


subject to the initial conditions z(0) = 1 = y(0), where a,b,c, and d are positive 
parameters 
e a: represents the natural growth rate of the pray in the absence of preda- 
tors; 
e b: represents the effect of predation of the pray; 
e c: represents the natural death rate of the predator in the absence of pray; 
e d: represents the efficiency and propagation rate of the predator in the 
present of pray. 
In numerical results we consider a = b = c= d = 1 and the implicit solution of the 
model is xy = exp(x + y — 2) 


5. CONCLUSION 
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